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ABSTRACT 



Context. To derive space velocities of stars along the line of sight from wavelength shifts in stellar spectra requires accounting for 
a number of second-order effects. For most stars, gravitational redshifts, convective blueshifts, and transverse stellar motion are the 
dominant contributors. 

Aims. We provide theoretical corrections for the net velocity shifts due to convection expected for the measurements from the Gaia 
■ Radial Velocity Spectrometer (RVS). 

Methods. We used a set of three-dimensional time-dependent simulations of stellar surface convection computed with C05BOLD to 
calculate spectra of late-type stars in the Gaia RVS range and to infer the net velocity offset that convective motions will induce in 
O ■ radial velocities derived by cross-correlation. 

Results. The net velocity shifts derived by cross-correlation depend both on the wavelength range and spectral resolution of the 
observations. Convective shifts for Gaia RVS observations are less than 0.1 km s for late-K-type stars, and they increase with stellar 
o3 ' mass, reaching about 0.3 km s 1 or more for early F-type dwarfs. This tendency is the result of an increase with effective temperature in 

both temperature and velocity fluctuations in the line-forming region. Our simulations also indicate that the net RVS convective shifts 
can be positive (i.e. redshifts) in some cases. Overall, the blueshifts weaken slightly with increasing surface gravity, and are enhanced 
at low metallicity. Gravitational redshifts amount up to 0.7 km s -1 and dominate convective blueshifts for dwarfs, but become much 
| weaker for giants. 

Key words. Stars: atmospheres - Stars: solar-type - Stars: late-type - Line: formation - convection - Techniques: radial velocities 

CO 

1. Introduction zone and the velocity field weakens with height. This causes 

' , . , , . . the net convective shift to be strongest for lines formed in deep 

CO; Stellar atmospheres are not in hydrostatic equilibrium, , progressively weakening and then disappearing for lines 
r-H , Convection, magnetic fields, rotation oscillations, and other formed M ^ ^ ^ hot here (Allende Priet0 & Garcfa 
> . phenomena induce time and spatially changing structure on the L6 1998; pierce & L Q 2000; RamfreZi Allende prieto 
surfaces of stars. In the case of late-type stars, surface convec- ^ L amrjer t 2008' Gray 2009) 
tion is responsible for the granulation pattern directly observed 
$-j . in high-resolution pictures of the Sun. Stellar spectra are often used to measure the velocities of 

Granulation broadens absorption line profiles in observed stars projected along the line of sight. The interpretation of 
solar spectra, making them asymmetric and introducing a net the measured line wavelength shifts in terms of Doppler shifts 
blueshift in their central wavelengths. Even though the typical purely associated with bulk stellar motion results in biased ra- 
granular velocities reach several kilometers per second, the net dial velocities owing to the presence of photospheric velocity 
blueshift observed in spatially-averaged spectra of solar photo- fields. Furthermore, motions in stellar atmospheres are not the 
spheric lines rarely exceeds 1 km s" 1 . In solar-type stars, the only mechanism that can bias radial velocity measurements, and 
photosphere is usually located above the convectively unstable a number of other effects, most notably gravitational redshifts 

~~ 7 (e.g., Pasquini et al. 2011), frustrate the direct derivation of 

The authors acknowledge the Texas Advanced Computing Center ^ u v r ■ , , , £ , £ , , , c 

rj,, j t . fc r _ . . j. r ... b TTri „ the line-of-sight component of the velocity of a star s center of 

(TACC) at The University ol Texas at Austin tor providing HPC ,. , r _ ■ r , n- ■ > 

resources that have contributed to the research results reported mass from the spectrum. CoiTecting for all these effects involves 

within this paper (www.tacc.utexas.edu). The model spectra dis- knowledge that is usually uncertain and, in most cases, model 

cussed in this paper are available in electronic form at the CDS dependent, and therefore the IAU recommends using the term 

via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 'radial- velocity measure' to refer to the velocity naively inferred 

http://cdsweb.u-strasbg.fr/cgi-bin/qcat ?J/A+ A/ from the spectral shifts under classical Newtonian mechanics, ig- 
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noring the corrections mentioned earlier (Lindegren & Dravins 
2003). 

The ESA mission Gaia is scheduled for launch in late 2013 
with the purpose of measuring proper motions, parallaxes, and 
spectrophotometry for some 10 9 stars down to V ~ 20 mag, and 
spectra for about 10 8 stars down to V ~ 17 mag (de Bruijne et 
al. 2009; Lindegren 2010). The typical uncertainty of the radial- 
velocity measures provided by Radial Velocity Spectrometer 
(RVS) onboard Gaia (Katz et al. 2004; Wilkinson et al. 2005) 
will typically be in the range 1-20 km s _1 , but taking the aver- 
age for large groups of stars will largely reduce the uncertainties, 
and the systematic errors discussed above can become dominant 
if neglected. More important, the RVS does not have any on- 
board calibration source to rely on for wavelength calibration, 
and systematic effects need to be considered, since wavelength 
calibration is based on a subset of the very same stellar spectra 
obtained by RVS (Katz et al. 201 1). 

As nature favors the formation of stars with low masses, we 
focus our attention on late-type stars, for which systematic errors 
affecting the translation from spectral shifts to bulk velocities are 
mainly gravitational redshifts, convective shifts, and the effect of 
the transversal motion on the relativistic Doppler effect. With the 
parallaxes measured by Gaia, through the comparison of colors 
and absolute brightness with stellar evolution theory, it will be 
possible to estimate masses and radii to calculate gravitational 
redshifts to better than 10 % (Allende Prieto & Lambert 1999), 
or some 50 m s _1 for a solar-like star. Accurate proper motions 
and parallaxes will also provide accurate transversal motions. 
With the exception of the Sun, it is very difficult to disentangle 
convective shifts from the bulk stellar motion, and therefore it 
becomes necessary to rely on models. 

We have shown that modern three-dimensional hydrodynam- 
ical models of the solar surface can accurately predict the con- 
vective shifts of individual weak and moderate-strength photo- 
spheric iron lines to within 70 m s _1 (Allende Prieto et al. 2009). 
The same models, however, predict net redshifts, which are not 
observed, for very strong iron lines (Allende Prieto et al. 2009; 
see also Asplund et al. 2000). Convective shifts need to be eval- 
uated for specific spectral windows used in the derivation of ra- 
dial velocities, ensuring that models are appropriate for those 
windows. 

In this paper, we focus our attention on the observations to 
be obtained by the Gaia RVS and predict convective wavelength 
shifts for stars of spectral types M through A (3500 < T e ff < 
7000 K) from a collection of hydrodynamical simulations pre- 
sented by Ludwig et al. (2009). As such corrections are specifi- 
cally derived for the Gaia RVS, and in particular for radial veloc- 
ity measures derived by cross-correlation of RVS spectra with 
template spectra calculated from classical (hydrostatic) model 
atmospheres, we will refer to them, indistinctively, as net con- 
vective shifts or effective convective shifts for RVS spectra. 

Section 2 describes the hydrodynamical simulations. Section 
3 is devoted to the spectral synthesis calculations, and §4 de- 
scribes the results. In §5 we comment briefly on the size of the 
gravitational shifts expected for the same stars previously mod- 
eled, and Section 6 presents a summary of our findings. 

2. Model atmospheres 

Simulations of stellar surface convection for 83 combinations of 
surface gravity, entropy flux at the bottom of the atmosphere, and 
chemical composition were computed using the code CO s BOLD 
(Freytag, Steffen & Dorch 2002; Wedemeyeret al. 2004; Freytag 
et al. 2012). The code solves the equations of hydrodynamics in 



a Cartesian grid with periodic horizontal boundary conditions 
and accounting for the effect of the radiation field on the energy 
balance. 

The models typically have a resolution of 140x 140x 150 (X- 
Y-Z) grid points, corresponding to physical sizes significantly 
larger than the typical size of granules, and range from a few 
Mm for dwarfs to thousands of Mm for giant stars. The radiation 
field applying a multi-group technique (Nordlund 1982) using 5 
groups for models of solar and 6 groups for models of subsolar 
metallicity. The number of groups was chosen by trial-and-error 
to ensure an acceptable representation of the radiative cooling or 
heating effects. This was done primarily by looking at solar-like 
stars. It turned out, however, that also giants, as well as hotter 
and cooler dwarfs could be treated in the same way obtaining 
similar accuracies. It was found that at subsolar metallicities the 
treatment of the radiative transfer in the continuum is more crit- 
ical, so that an extra group was introduced. 

We believe that the typically employed five or six groups are 
sufficient for obtaining accurate line shifts, but this is difficult to 
demonstrate in general since only for few cases are models with 
more groups available. However, the line shifts hinge very much 
on the velocity field which is driven in deeper layers, where the 
wavelength dependence of the radiative transfer becomes of sec- 
ondary importance. In the few models employing 12 groups that 
are available, changes in the level of the temperature fluctuations 
are not large. We think that rather spatial resolution and sampling 
of the time series dominate the uncertainties in the line shifts. 

The equation of state takes into account H and He ioniza- 
tion and H2 formation, while the opacities are calculated with 
the Uppsala package (Gustafsson et al. 2008). The physical 
size of the computational box is the same for models with the 
same surface temperature and gravity. The upper boundary of 
the models is open, and extended beyond a Rosseland optical 
depth of 10 -6 , while the lower boundary is well into optically 
thick, convectively unstable layers. Each simulation was run un- 
til fully relaxed. Chemical abundances are from Grevesse & 
Sauval (1998), with the exception of CNO, which are updated 
following Asplund (2005). More details are provided by Ludwig 
etal. (2009). 

For each simulation, a set of uncorrelated snapshots was se- 
lected for spectral synthesis. The number of snapshots and the 
total time span of the simulations varied from star to star, ac- 
cording to the granular turn-over time scale for each model. The 
snapshots were hand chosen in order to be representative of the 
whole simulated time series, avoiding statistical anomalies. We 
found from tests that about 20 uncorrelated snapshots are enough 
to represent a model, but in some cases a smaller number was 
adequate. 

In order to determine the net convective shift predicted for 
each 3D model, we calculated a ID model atmosphere with the 
same effective temperature, surface gravity and metallicity. The 
effective temperature for each simulation was determined from 
the temporally and spatially averaged emergent bolometric flux. 
The ID models were calculated with the LHD package (Caffau 
& Ludwig 2007), and share the basic ingredients (most impor- 
tantly the opacities and equation of state) with the 3D simula- 
tions. A second comparison set of ID models was derived by 
linear interpolation of the ODFNEW model grid by Castelli & 
Kuruc20 (2004), based on the solar abundances of Grevesse & 
Sauval (1998). 
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Fig. 1. Observed (black line) and synthetic (ID in red, 3D in blue) spectra for Arcturus and the Sun - all data smoothed with a 
Gaussian kernel to have a resolving power of about 11,500 as appropriate for the Gaia RVS. Both the normalized spectra and the 
residuals are shown for each star. A microturbulence of 1.5 km s was adopted for the ID models, which were calculated with the 
LHD package. 



3. Spectral synthesis 

3.1. 3D spectral synthesis 

We computed synthetic spectra covering the Gaia RVS wave- 
length range (847-874 nm) with wavelength steps equivalent to 
< 0.1 km s _1 (or a resolving power R = A/6A ~ 10 6 ). The cal- 
culations were carried out with the 3D synthesis code ASSeT 
(Koesterke 2009; Koesterke, Allende Prieto & Lambert 2008). 
The reference solar abundances are from Asplund, Grevesse & 
Sauval (2005), and while temperature, density and velocity fields 
were taken from the model atmospheres, the electron density 
was recalculated with an equation of state including the first 99 
elements in the periodic table and 338 molecules (Tsuji 1973, 
with partition functions from Irwin 1981). 

The number of snapshots used varied between 8 and 2 1 , but 
in the vast majority of cases was about 20, and the total time 
covered ranged between 0.6 and 19500 hours, spanning at least 
10 times the typical lifetimes of convective cells. The number of 
pixels considered for each snapshot was reduced by taking one 
out of three points in the X- and Y-directions, and at least the 
uppermost 90 vertical points. These choices are based on previ- 
ous experience and limited testing with the set of models used 
in this paper. These tests indicate that dropping spatial points 
would lead only to differences at the level of about 0.2% for a 



typical strong Fe I line. For the continuum and weak lines the 
differences tend to be much smaller than that. 

The adopted solar abundances are very similar to those used 
in the 3D simulations. The small existing inconsistencies are ex- 
pected to have a negligible impact in the predicted convective 
shifts. Continuous absorption from H and H is included. Line 
absorption is included in detail from the atomic and molecu- 
lar files compiled by Kurucz. Opacities are precomputed in a 
temperature-density grid bracketing all the snapshots with steps 
of 250 K and 0.25 dex, and subsequently derived for the grid 
points using cubic Bezier interpolation. 

For each individual snapshot a spectrum is derived in 2 steps. 
First a background radiation field is calculated. This is used for 
the scattering term in the subsequent spectrum synthesis. Note 
that actually 2 spectra are calculated for each snapshot. The sec- 
ond one, which is derived from opacity data without spectral 
lines, is used for continuum normalization. 

The radiative transfer calculations account for atomic hydro- 
gen Rayleigh and electron scattering. The mean background ra- 
diation field is calculated from 24 angles (8 azimuthal with lin- 
ear weights, and 3 zenith per half-sphere with Gaussian weights) 
using a short characteristics scheme, while the emergent flux is 
integrated for 2 1 angles (vertical ray plus 3 zenithal rays and 8 
azimuthal angles again with linear weights, except that only 4 
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zenithal angles are considered for the shallowest rays). The ab- 
scissae and weights for each case can be found in Abramowitz & 
Stegun (1972). This setup evenly covers a sphere with 8 angles 
for the azimuth (0-360 deg) and 3-4 angles for the inclination 
(0-180 deg). The integration along a ray was performed from the 
top layer down to optical depths of at least 20. The final emer- 
gent flux is an average of the flux for all the snapshots considered 
for each simulation. 

Rotation was not accounted for in this step. It is well-known 
that rotation affects the observed line asymmetries (Gray 1986), 
but low rotational velocities should only have a modest effect 
on the net shifts derived by cross-correlation for the entire Gaia 
RVS spectral window. In part, this is related to the spectral res- 
olution of the RVS instrument, corresponding to 26 km s _1 (de 
Bruijne 2012). 

The calculations were performed on the Longhom cluster at 
the Texas Advanced Computing Center (TACC). The intensities 
for about 150 million frequencies were derived across all snap- 
shots. The total computational effort is roughly equivalent to 3 
months on a modern workstation with 4 cores. 

Fig. Q] illustrates the spectra in our ID and 3D grids that 
are closest to the atmospheric parameters for Arcturus and the 
Sun. This comparison is done after smoothing the data to a 
FWHM resolving power of R - 11,500 as appropriate for 
the RVS; we adopted a microturbulent velocity of 1.5 km s _1 
for the ID calculations. The agreement between models and 
observations is significantly better in the solar case (observed 
spectrum from Kurucz et al. 1984) than for Arcturus (Hinkle 
et al. 2000). Overall, predicted absorption features are slightly 
stronger in ID than in 3D. The agreement is reasonable, consid- 
ering that the models are not tailored to these stars. Arcturus has 
an r e ff ~ 4300 K, logg s 1,7 (with g in cm s~ 2 ) and a metal - 
licity about [Fe/H] ^ -0.5 (see, e.g., Ramirez & Allende Prieto 
201 1), and the Sun has an T eS = 5777 K, logg 4.437 and, by 
definition, [Fe/H]= (see, e.g. Stix 2004). The model compared 
here with Arcturus (d3t40gl5mml0n01) has an r eff of 4040 K, 
\ogg =1.5 and a metallicity of [Fe/H]= —1.0, and that com- 
pared with the Sun (d3t59g45mm00n01) has an T eS of 5865 K, 
log g = 4.5 and solar metallicity. 

3.2. 1 D spectral synthesis 

The ID counterpart radiative transfer calculations were per- 
formed with the ID version of ASSeT, and therefore the results 
can be directly compared between the ID and 3D models. In 
particular, chemical compositions, opacities, and radiative trans- 
fer calculation for the emergent intensities are shared by the ID 
and 3D branches. The only differences between the ID and 3D 
calculations are that the opacities are calculated exactly for the 
grid points in ID (every depth in the model atmospheres), while 
they are interpolated in 3D (from a precomputed grid in tem- 
perature and density, as described in 33. Il l, and the radiative 
transfer solver for the calculation of the mean intensities in the 
grid, which is based on Feautrier's method (Feautrier 1964) in 
ID and an integral method in 3D. Our tests nonetheless indicate 
that these differences have a negligible effect on the continuum- 
normalized fluxes. The difference between the Feautrier and the 
direct method in ID is about 10~ 5 . The effect of the interpolation 
of the opacities in 3D is estimated to be on the order of 10 4 . The 
different choices for the opacity calculation and the radiation 
transfer method are made to accommodate the very expensive 
calculations in 3D, which are not feasible without interpolating 
the opacities and benefit from the faster 1 st order method for the 
integration of the intensities. 
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Fig. 2. The spectra for a single snapshot of model 
d3t59g45mm00n01 (black) and its ID counterpart (black), 
with near solar atmospheric parameters, are shown in the 
left-hand panels for two different values of the resolving power 
(R). The right-hand panels illustrate how the corresponding 
shifts (ID relative to 3D) are obtained by fitting the peak of the 
cross-correlation function with a parabola. 

4. Effective convective shifts for the Gaia RVS 

We divided the spectra by the pseudo continua, estimated by it- 
eratively fitting a low order polynomial to the upper envelope of 
each spectrum. We then calculated the net offset between the 3D 
and ID spectra to derive the overall convective shift predicted by 
the hydrodynamical simulations by cross-correlation. 

The cross-correlation of two arrays (o spectra) T and S of 
equal and even number of elements N is defined as a new array 
C with 



C = J] T k S 



(1) 



where i runs from 1 to N. If the spectrum T is identical to S, 
but shifted by an an integer number of pixels p, the maximum 
value in the array C will correspond to its element i = p + j. 
Cross-correlation can be similarly used to measure shifts that 
correspond to non-integer numbers. In this case, finding the lo- 
cation of the maximum value of the cross-correlation function 
can be performed with a vast choice of algorithms. 

We calculated the cross-correlation of ID and 3D model 
spectra using the software by Allende Prieto (2007), which al- 
lows fitting the peak of the cross-correlation function with a 
parabola, a cubic polynomial, or a Gaussian. Based on tests 
adding noise to simulated spectra, we chose to use parabolic fit- 
tings to the central 3 points around the maximum of the cross- 
correlation function. 

The net convective shift derived by cross-correlation not only 
depends on the chosen spectral range, but also on the resolv- 
ing power of the spectrograph. A degraded resolution will in- 
troduce, in addition to random noise, a systematic error in the 
derived velocity shifts. This is illustrated in Fig. [2] for a single 
3D snapshot and its ID counterpart at two different values of 
the resolving power: resolution does change the apparent con- 
vective shift. Fig.|3]shows the derived shifts (ID relative to 3D) 
for model d3t59g45mm00n01 (see Table 1 for reference), which 
was compared with the solar spectrum in Fig. Q] as a function 
of the resolving power. At the resolving power of the Gaia RVS 
(R =11,500), the inferred convective blueshift is about 0.1 km 
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FWHM resolving power (R) 

Fig. 3. Derived convective shift for model d3t59g45mm00n01, 
with near solar atmospheric parameters, as a function of the re- 
solving power. These shifts are obtained by cross-correlation us- 
ing the Gaia RVS spectral window, and fitting the peak of the 
cross-correlation function with a parabola. The dashed vertical 
line marks the RVS resolving power. 

s _1 , while that measured when the spectral lines are fully re- 
solved (R > 200,000) approaches 0.24 km s . Our results for 
the simulations are included in Table 1 . All the values presented 
in the table, and discussed below correspond to the RVS resolu- 
tion. 

In addition to calculating the net convective shifts for the av- 
erage spectra from the selected snapshots for each simulation, 
we also calculated the shifts for each snapshot and averaged 
them out. These correspond to the 'mean' shifts in Table 1. The 
variance (cr 2 ) among snapshots was used to estimate the intrinsic 
uncertainties as cr/ V«> with n the number of snapshots consid- 
ered, and these are also included in the table. We believe that the 
uncertainties calculated in this way are conservative since the 
subsample of snapshots selected for spectral synthesis is chosen 
to represent the properties of the complete run (e.g., in effec- 
tive temperature and mean velocities). Moreover, the separation 
in time among the snapshots is large enough that they can be 
considered as independent despite being selected from one time 
series. The agreement between the shifts of the average spec- 
tra and the average shifts from individual snapshots is excellent, 
with an average of 0.0004 km s _1 and a standard deviation of 
0.001 km s" 1 . 

Fig. [4] shows the convective shifts for all models. Note that 
there are models for four different metallicities, color-coded, and 
that the resulting effective temperatures of the simulations are 
not always exactly the same for any combination of surface grav- 
ity and metallicity, as this is not an input parameter - the entropy 
of the material entering the simulations box through the open 
bottom boundary is a chosen input instead. 

The changes of the net convective shifts with the atmospheric 
parameters are perhaps most obvious in Fig. |5] The net shifts 
show a mild correlation with the effective temperature of the 
models, in the sense that warmer stars tend to exhibit larger 
blueshifts. As found in previous studies (Dravins & Nordlund 
1990a,b; Nordlund & Dravins 1990), the net energy flux travers- 
ing the atmosphere is the most important parameter. Surface 
gravity has a small effect on the convective shifts predicted by 
our simulations - overall the shifts intensify at lower gravities 
-, and so has metallicity, with lower net blueshifts at higher 



metallicities. Fig. |4] shows the results obtained when the LHD 
ID models, which are fully consistent with the 3D simulations 
in terms of opacities and the equation of state, are used as ref- 
erence. However, very similar results are found when Kurucz 
models are used as a reference (see below). 

We adopted a micro-turbulence of 2.0 km s _1 , and a macro- 
turbulence of 1.5 km s _1 for the reference ID calculations, 
but our experiments indicate that these parameters would only 
lead to a modest impact on the spectral shifts determined by 
cross-correlation. For instance, completely neglecting macro- 
turbulence leads to a mean correction in the predicted convective 
blueshifts of just 0.01 (cr = 0.032) km s _1 . Adopting a micro- 
turbulence of 1.5 km s _1 instead of 2.0 km s _1 will, on average, 
reduce the convective blueshifts by about 0.007 (cr = 0.011) 
km s _1 . Finally, adopting Kurucz model atmospheres instead of 
LHD models as our ID reference will, on average, enhance the 
blueshifts by about 0.01 (cr = 0.03) km s" 1 . 

The smallest (in absolute terms) convective shifts in the grid 
are for the coolest dwarfs of low metal content, and amount 
barely to < 0.1 km s _1 . The most vigorous net velocities 
are found for low metallicity F-type subgiants, practically the 
warmest stars considered in the grid, where they exceed 0.3 km 
s- 1 . 

There are a few simulations (mostly with [Fe/H]= -1) for 
which redshifts are predicted. While perhaps not very evident 
in Fig. [5] the drop is clearly apparent in a global analytical fit 
of the convective shifts of the 3D simulations (see AppendixfAl. 
These redshifts are a clear prediction of the simulations. They are 
caused by strong redshifts of the cores of the Ca II triplet lines, 
but they dominate the overall cross-correlation signal. The cores 
get redshifted when the lines are strong and the flow exhibits a 
pronounced "reverse granulation" pattern. Still, the overall shift 
of the Ca II triplet line cores is a superposition of redshifted and 
blueshifted contributions and therefore it is difficult to predict a 
priori which ones are dominating. 

We tried to interpret our findings in terms of the temperature 
and velocity fluctuations present in the 3D simulations. We aver- 
aged the temperature and vertical velocity component horizon- 
tally and temporally on surfaces of constant Rosseland optical 
depth. We took values at Rosseland depth t = 2/3 as represen- 
tative for the conditions prevailing in the line forming regions - 
at least for weak lines. These are shown in the two bottom pan- 
els of Fig. |5] where a mild correlation between temperature and 
velocity fluctuations is apparent. 

Figure [6] shows how the temperature and velocity fluctua- 
tions depend on effective temperature, surface gravity and metal- 
licity. Interestingly, we find that the velocity fluctuations depend 
on metallicity for T e ff < 5000 K but they do not for warmer 
stars. Only close to the main-sequence a systematic dependence 
appears below 5000 K with higher metallicity implying greater 
fluctuations. For main-sequence stars, r e jf = 5000 K also sep- 
arates regions of different correlation between metallicity and 
temperature fluctuations. Above 5000 K, higher metallicity cor- 
responds to lower temperature fluctuations but this is reversed 
below 5000 K. While Fig. [6] certainly motivates the increase 
in convective shifts towards higher effective temperatures and 
lower gravity, it cannot explain all trends of the convective shifts. 
What is indeed missing is information on the mutual correla- 
tion between temperature and velocity fluctuations, as well as 
on changes of the line formation conditions with atmospheric pa- 
rameters - as discussed further below. In particular, Fig. |6]does 
not explain the fast drop of the predicted convective shifts near 
the highest effective temperatures in our grid. 
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Fig. 4. Effective convective shifts predicted by the 3D hydrodynamical models for Gaia RVS observations. 



Nevertheless, our results resemble earlier high-resolution 
studies of spectral line shapes and convective shifts, but the net 
velocity shifts predicted for the RVS observations are smaller 
than for individual, resolved, spectral lines (see, e.g., Bigot & 
Thevenin 2006, 2008; Chiavassa et al. 2011). As expected, the 
net shift measured by cross-correlation is a weighted average of 
the lineshifts from all spectral features. At high metallicity, many 
lines contribute, with strong lines giving very small convective 
shifts, and weaker features contributing larger shifts. For exam- 
ple, for a solar-like atmosphere, if we apply cross-correlation 
by isolating windows around the Ca II lines, we find very small 
convective shifts on the order of 0.030 km s _1 , while limiting the 
spectral window to the stretch between the Ca II lines (857-863 
nm), the net shift is about -0.18 km s _1 , double than what we 
obtain for the whole RVS spectral range. 

For a solar 3D model not included in the grid discussed here, 
we found in a previous paper (Allende Prieto et al. 2009) a con- 
vective shift of -0.26 km s , in excellent agreement with the 
value directly derived from the solar flux atlas of Kurucz et al. 
(1984). These two values correspond to observations at much 
higher spectral resolution, and can be compared with those we 
find (see Fig. 2) for the model in the grid with parameters closest 
to solar, d3t59g45mm00n01, namely -0.24 km s when the 
resolving power is R > 200, 000 but note that changes to -0.09 
km s _1 atR = 11,500. 

From a theoretical standpoint, the main effect of a reduced 
gravity is a lower surface pressure, as well as an increase in both 
the vertical and horizontal scales. At a given effective temper- 
ature the flow is constrained to transport the same energy flux 
so that a reduced pressure and density (since the temperature is 
almost fixed) has to be compensated for by an increase in the ve- 



locity and/or the temperature fluctuations (Dravins & Nordlund 
1990b; Dravins et al. 1993; Collet, Asplund & Trampedach 
2007; Gray 2009). As a result, convective blueshifts should 
strenghten - as long as the surface area fractions covered by up- 
flows and down-flows do not change. As illustrated in Fig. [6] our 
detailed simulations roughly support this simple picture. 

At lower metallicity we expect the reduced continuum opac- 
ity to make visible deeper atmospheric layers with more vig- 
orous convection (Allende Prieto et al. 1999), and the simula- 
tions indicate that the convective shifts of individual lines be- 
come stronger, but the span of the line bisectors weakens dra- 
matically with metallicity. This is illustrated for the strong Fe I 
line at 868.86 nm in solar-like (T eff =* 5900 K and logg = 4.5) 
stars in Fig. [7] In this figure, we show line profiles (top panel) 
and the corresponding bisectors (lower panel), with the solid 
lines for the calculations based on 3D models and the dashed 
lines for ID models. This analysis is performed at full resolu- 
tion, since line bisectors are erased at the RVS resolution, using 
a micro-turbulence of 2 km s _1 and no macro-turbulence. As the 
feature weakens, the range of atmospheric heights sampled from 
the wings to the core is reduced, which is a likely explanation 
for the smaller bisector spans. For this line, we also find that the 
predicted convective shift is -0.00, -0.10, -0.12, and -0.13 km 
s _1 at [Fe/H]= 0,-1,-2, and -3, respectively. Note that these 
bisectors are calculated using the standard definition, i.e. sub- 
tracting the velocity in the blue wing from that at the same nor- 
malized flux level in the red wing, and thus the net line shifts 
are not included and they approach zero velocity in the line core 
(Gray 1992). In addition to the characteristics of surface convec- 
tion changing as a function of metallicity, the strengthening of 
the line blueshifts with decreasing metallicity is related to the 



6 



C. Allende Prieto et al.: Convective line shifts for the Gaia RVS from the CIFIST 3D model atmosphere grid 



X 
CD 



6500 
, 6000 
• 5500 
5000 
4500 
4000 
4.5 
4.0 

r 3.5 

3 3.0 
cp 2.5 

CP 

° 2.0 

1 .5 
0.0 

-0.5 

-1 .0 

-1 .5 

-2.0 
-2.5 
-3.0 
0.06 
0.05 
0.04 
0.03 
0.02 

0.01 

0.00 

P 4 
w 

E 3 
X 2 
1 





E 



-0.4 




• • ••• 



• , — f mm f ■ f» 



• • 

• * 5 • 



-0.2 0.0 

Convective shift (km 



0.2 



0.4 



Fig. 5. Three uppermost panels: Predicted convective shifts and atmospheric parameters of the 3D hydrodynamical models for Gaia 
RVS observations. Two lowermost panels: Convective shifts versus relative horizontal temperature fluctuations and the rms of the 
vertical velocity at Rosseland optical depth r — 2/3. In all panels data points are color coded by metallicity as in Fig. [4] 



line formation shifting deeper into the photosphere as the line 
weakens. This is consistent with our findings for the overall net 
blueshift for Gaia RVS observations. 

The filled circles in the upper panel of Fig. |7]show the renor- 
malized observations for the Sun of Kurucz et al. (1984). A sig- 
nificant improvement in the agreement between theory and ob- 
servation is obvious for the line core when going from ID to 3D 
models (dashed and solid black lines, respectively). A similar 
effect is found using a Kurucz model, instead of its LHD coun- 
terpart, for the ID spectrum. It is also noticeable the extreme 
discrepancies between the ID and 3D cases at [Fe/H]= -2, as 
it is the large reduction in strength for the line calculated in 3D 
between [Fe/H]= -2 and -3. An exhaustive check of these pre- 



dictions against high quality observations of metal-poor stars is 
very much needed. 

Our calculations are likely least reliable for cool giants, and 
in fact the predicted convective shifts for stars approaching 4000 
K are unlikely to be realistic. Gray, Carney & Yong (2008) ex- 
amined bisectors for a number of red giant stars with low metal- 
licity, finding that those with temperatures of about 4100 K or 
less tend to exhibit reversed bisectors (with an 'inverse C shape, 
rather than with a 'C'-like shape seen in solar-type stars). More 
recent results for the supergiant y Cyg confirm the same trend 
(Gray 2010). These stars also tend to exhibit radial velocity 'jit- 
ter' (see, e.g., Carney et al. 2008, and the theoretical predictions 
by Chiavassa et al. 201 1). 
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Fig. 6. Relative horizontal temperature fluctuation (top panel) 
and absolute vertical velocity fluctuations (bottom panel) at 
Rosseland optical depth t = 2/3 in the 3D models. The sur- 
face gravity is color coded, metallicity is depicted by different 
plot symbols. 



We have identified an expression that can reproduce the net 
convective blueshifts theoretically predicted for the RVS obser- 
vations in Table 1 with an rms of 0.021 km s _1 ; and a maximum 
deviation of 0.074 km s _1 . We have coded this expression into 
an IDL function provided in the appendix. 



5. Gravitational redshifts 

Photons escaping from a stellar photosphere to infinity are red- 
shifted as a result of climbing up the gravitational field by 
GcM/R, where G is the Gravitational constant, c the speed of 
light, M the stellar mass, and R the stellar radius. Using solar 
units for the stellar mass and radius, the gravitational redshift 
amounts to 0.6365 M/R km s _1 , and for light intercepted at 1 
AU 0.6335 M/R km s" 1 (Lindegren & Dravins 2003). The same 
effect acts when photons enter in the Earth's gravitational field, 
which reduces the redshift by about 0.2 m s _1 , but this small 
correction can be ignored for most purposes. With the parallaxes 
and spectrophotometry provided by Gaia, it should be possible 
to obtain fairly accurate customized predictions for the gravita- 
tional redshifts for individual Gaia targets. 

Gravitational shifts for late-type dwarfs amount to about 0.7- 
0.8 km s at \ogg ~ 4.5, but they dramatically decrease with 
surface gravity down to 0.02-0.03 km s _1 for the giants in our 
grid (logg -1.5). 
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Fig. 7. Line profiles and velocity bisectors for the Fe I transition 
at 868.6 nm. The filled circles in the upper panel correspond to 
solar observations. Solid lines are for 3D simulations and bro- 
ken lines for ID models. Note that the ID calculation at solar- 
metallicity for this line (broken solid line) shows a non-zero ve- 
locity bisector, due to overlapping CN transitions included in all 
our calculations that distort the otherwise perfectly symmetric 
ID profiles. 



6. Discussion and Summary 

The radial velocities inferred from the wavelength shifts mea- 
sured in stellar spectra are systematically offset from the true 
space velocities of stars projected along the line of sight due to 
a number of effects. Quantifying the impact of such effects is 
model dependent, and therefore the IAU has introduced the term 
'radial- velocity measure' (Lindegren & Dravins 2003) to refer to 
the radial velocity naively inferred from the relative wavelength 
shifts of lines c6A/A, and encourages observers to keep measure- 
ments disentangled from subsequent correctionsQ. Nevertheless, 
for most stars, two relatively clean corrections are dominant: 
photospheric gravitational redshifts and convective blueshifts. In 
the case of stars moving at high speed in the direction perpen- 
dicular to the observer, the relativistic contribution to the Dopper 
shift needs to be considered at a precision level of hundred of 



2 We note that the 'radial-velocity measure' does not have a unique 
definition, since the value derived depends on the spectral features and 
the actual method used. Therefore, such measures must be accompanied 
by a detailed description of the procedure used to derive it. 
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meters per second; the spectrum of a star with a transverse ve- 
locity of 300 km s _1 and no radial velocity will be redshifted as 
much as that of a star with a receding radial velocity of 0. 15 km 
s" 1 . 

We have used three-dimensional radiative-hydrodynamical 
simulations of surface convection in late-type stars to evaluate 
the impact of convective shifts on the spectra from the Radial 
Velocity Spectrometer onboard Gaia. The net velocity shifts de- 
rived by cross-correlation depend both on the wavelength range 
and the spectral resolution of the observations. Our models pre- 
dict very small net convective shifts for the RVS observations of 
K-type stars, which increase with effective temperature to reach 
bout 0.3 km s _1 for F-type dwarfs. The predicted shifts depend 
as well on metallicity and surface gravity. 

In the case of the Gaia mission, the corrections discussed 
above, which are < 1 km s are small in comparison with the 
typical uncertainties in the space velocities determined for indi- 
vidual stars, but can cause systematic errors when average values 
are derived for ensembles of stars. Even more important, since 
the Gaia RVS is a self-calibrated instrument, in the sense that the 
wavelength solution is derived from stellar observations, system- 
atic offsets from star to star such as those related to convective 
shifts, need to be accounted for in the data reduction. In this 
paper, we used a grid of 3D hydrodynamical simulations to es- 
timate convective blueshifts for stars with spectral types F to M, 
at different metallicities and evolutionary stages. We note that 
from data provided by Gaia itself, it will be possible to constrain 
the stellar masses and radii for individual stars, and therefore the 
gravitational redshifts. 

Gravitational redshifts dominate convective blueshifts pre- 
dicted for the Gaia window for dwarfs of all temperatures, in 
particular for later spectral types (K and M), where the convec- 
tive shifts are weakest. Gravitational redshifts become weaker 
for stars with lower gravities. The opposite tendency, if less pro- 
nounced, is expected for the convective blueshifts, in such a way 
that gravitational and convective shifts nearly cancel each other 
for stars with intermediate gravities. 

When adding together convective and gravitational shifts, 
we expect the largest velocity offset for the coolest dwarfs in 
our study, with effective temperatures around 4000 K, where the 
net effect, dominated by gravitational redshifts, amounts up to 
nearly 0.6 km s _1 . By using the corrections provided here and 
custom-made estimates of the gravitational shifts, these sources 
of systematic error should be reduced significantly to a level 
ranging from about 0.1 km s" 1 for F-type dwarfs, to a few tens 
of meters per second for late-K-type stars. 

In addition to calculated convective shifts, we provide FITS 
files including the ID and 3D model fluxes we have computed. 
These may turn useful to build radial velocity templates for 
derivation of radial velocities from observations. 

We end with a word of caution. The theoretical calculations 
presented here only have been checked against a limited set of 
observations, mainly for the Sun. An exhaustive check against 
high-quality observed spectra for warmer and cooler tempera- 
tures, especially at very low metallicity, is still needed before 
one can ensure that the predicted shifts are completely trustwor- 
thy. 
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Appendix A: An empirical fitting function to 
describe the blueshifts as a function of 
atmospheric parmeters 

To derive from the tabulated data an easily applicable function 
we follow the recursive fitting procedure outlined by Sbordone et 
al. (2010): we searched for an analytical function for the convec- 
tive blueshift with three independent parameters (log T e ff , log g 
and [Fe/H]), written as y — y(xo, X\, Xi\ A). The function and the 
vector of free parameters A = {Aq,A\, . . .) that give the best fit 
to the tabulated data is to be found. 

The function-find routine ("fufi", written in IDL) starts with 
the most simple "function", a constant Aq - the weighted average 
of the tabulated blueshift data. At each recursive iteration we let 
the routine replace each parameter A, of a candidate fitting func- 
tion with a polynomial A,- — > Aq +A\Xq +A2X1 +A3X2 or, alterna- 
tively, an exponential A, — > Aq+Ai exp(A2 -1-A3.ro +A5X2) 
to reach the next level of complexity. For each of these func- 
tions, the optimum parameter set is then searched, using the pa- 
rameters from the previous level as starting point. We use the 
inverse of the rms fluctuations of the convective shifts between 
the snapshots for a given stellar parameter combination (or 0.01 
km s-1 when the fluctuations are below this value) as fitting 
weights. The function with the fitting-parameter set that provides 
the smallest overall error (among thousands of candidates) is fi- 
nally written to a file as a Fortran or IDL function. The function 
we supply has been slightly edited to improve readability. We 
emphasize that the functional form and fitting parameters have 
no physical meaning whatsoever. 

We varied some control parameters (the list of candidate 
terms and the recursion depth), and compared errors and the 
overall functional dependence of the results. Fitting functions 
composed exclusively of simple polynomial terms already give 
decent fits (and would require a much simpler procedure to de- 
rive than the one outlined above). Still, fits including the ex- 
ponential term are superior and behave much better and closer 
to our expectations at - and even slightly beyond - the bor- 
ders of the set of stellar parameters for which model results are 
available. The overall rms scatter between the direct measure- 
ments from the simulations and the fitting formula is just 0.021 
km s _1 , with a maximum deviation of 0.074 km s _1 , We give 
an implementation of the fitting function in IDL below. It in- 
cludes a check that the input parameters are within the range 
covered by the 3D models. A few of the free fitting parameters 
are zero, but are left in the routine to allow for a more systematic 
writing of the final expression. Figs. IA.1I IA.2I IA.3I [Aj4] illus- 
trate the convective shifts predicted by the fitting equation for 
[Fe/H]= 0, -1,-2, and -3, respectively. 



; (<0 for blueshifts) 

; --- Thu Oct 20 20:09:5(8 2011 --- B.F. --- 
; check that 3 parameters are input 
if n_pararas() It 3 then begin 
print, ' 

return , - 10000000000 . dO 
endif 

; check that we're within limits 
if (min(Teff) It 3790.00 or $ 
raax(Teff) gt 6730.00 or $ 
rain(FeH) It -3.0 or $ 
raax(FeH) gt 0. or $ 
max(logg) gt 4.5 or $ 

min(logg-(teff*9. 184e-4-2 . 482)) It 0.0) $ 
then begin 

print, '% ConvShift: Params. out of range:' 
print, '% 3790 <= Teff <= 6730. K' 
print, '% 9. 184e-4*Teff-2 .482<=logg<= 4.5' 
print, '% -3.0 <= [Fe/H] <= 0.0 dex' 
return , - 10000000000 . d0 
endif 

logTeff=alogl0(Teff*ld-3) 
M=-FeH 

A=dindgen(24) 

A(00)= 3.9921954E-01 & A(01)=-7 . 1438992E-01 
A(02)= 2.0028012E-02 & A(03)=-5 . 2 146580E-02 
A(04)=-4. 1100047E-06 & A(05)= . 0000000E+00 
A(06)= 2.2334658E+01 & A(07)=- 1 . 70702 15E+00 
A(08)=-1.7445524E-01 & A(09)= 1 . 4332379E-02 
A(10)= 0.0000000E+00 & A( 1 1)=- 1 . 3239030E+01 
A(12)= 8.2077384E-01 & A(13)=-9.6311294E-©2 
A(14)=-1.7150490E+03 & A(15)= . 0000000E+00 
A(16)=-9.3405848E+00 & A(17)= 7 . 7603954E-01 
A(18)= 1.3519228E+00 & A(19)=-5 .0111694E+00 
A(20)= 0.0000000E+00 & A(2 l)=-8 . 2077414E-01 
A(22)= 1.0457402E-01 & A(23)= 4 . 2118204E-01 
Shift=A(00)+A(01)*logTeff+(A(02) + $ 
(A(04)+A(09)*exp(A(10)+(A(ll)+ $ 
A(14)*exp(A(15)+ $ 
(A(16)+A(19)*exp(A(20)+ $ 
A(21)*logTeff+ $ 

A(22)*logg+ A(23)*M))*logTeff+ $ 

A(17)*logg+A(18)*M))*logTeff+ $ 

A(12)*logg+A(13)*M))*exp(A(05)+ $ 

A(06)"logTeff+A(07)-logg+A(08)*M))*logg+A(03)*M 

Shifts- Shi ft 

return, Shift 

end 



function ConvShift, Teff, logg, FeH 

; IN: Teff -- effective temperature (K) 

; logg -- surface gravity (g in cm/s2) 

; FeH -- [Fe/H] metallicity relative to solar 

; OUT: Convective shift (m/s) 
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Fig.A.l. Convective shifts for solar metallicity simulations (in Fig. A.4. Convective shifts, [M/H]=-3. Similar to Fig. IA.ll 
units of m s _1 ). The contours join the locations with the same 
convective blueshift, as predicted by our analytical fitting func- 
tion. Red/green color is used in the regions where the analytical 
function under/over predicts the measurements performed on the 
simulations. 
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Fig. A.2. Convective shifts, [M/H]=-l. Similar to Fig. IA.ll 
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Fig. A.3. Convective shifts, [M/H]=-2. Similar to Fig. IA.ll 
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Table A.l. Predicted convective line shifts. The fifth column gives the convective 
shifts measured in the average spectrum from all snapshots considered, and the 
sixth column gives the average of the shifts measured in individual snapshots. 



Simulation ID T e jf logg [Fe/H] Convective shift Mean Error 

(K) (cms- 2 ) (dex) (kms^ 1 ) (kms^ 1 ) (kms^ 1 ) 



d3t38g40mm00n2 


3813 


4.00 


-0 


.00 


-0 


,043 


-0 


.043 





,004 


d3t38g40mml0n3 


3790 


4.00 


-1 


,00 


-0 


,019 


-0 


.019 





.001 


d3t40gl5mm00n02 


4018 


1.50 


-0 


.00 





,000 





,005 





.054 


d3t40gl5mml0n01 


4040 


1.50 


-1 


.00 





,123 





,124 
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